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Modeling excitable systems: Reentrant tachycardia 
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Excitable membranes are an important type of nonlinear dynamical system and their study can 
be used to provide a connection between physical and biological circuits. We discuss two models of 
excitable membranes important in cardiac and neural tissues. One model is based on the Fitzhugh- 
Nagumo equations and the other is based on a three-transistor excitable circuit. We construct a 
circuit that simulates reentrant tachycardia and its treatment by surgical ablation. This project 
is appropriate for advanced undergraduates as a laboratory capstone project, or as a senior thesis 
or honors project, and can also be a collaborative project, with one student responsible for the 
computational predictions and another for the circuit construction and measurements. 



I. INTRODUCTION 

Many physics majors are unaware of the important role 
that physicists play in medical research. Typical course- 
work for a physics major is removed from real- world ap- 
plications, especially compared to what they perceive for 
their pre-medical peers. As a result, some students ques- 
tion whether majoring in physics will allow them to pur- 
sue a useful and satisfying career. These students will 
benefit from participation in projects that make connec- 
tions between physical concepts and real-world medical 
problems. In this paper we focus on the modeling of 
reentrant tachycardias as a means to introduce physics 
majors to the physical underpinnings of excitable tissues. 

The voltage changes that occur across cell membranes 
during action potentials (the signal pulses) in cardiac and 
neural tissues constitute an important type of excitable 
nonlinear behavior. Excitable systems typically display 
a threshold-triggered fast-response positive feedback fol- 
lowed by a slow-response inhibition that causes the sys- 
tem to return to its resting values. The recovery includes 
a phase known as the refractory period during which ad- 
ditional excitation is either more difficult or impossible 
to elicit. Excitation is manifested by the temporal pulse 
shape of a system variable such as the membrane po- 
tential during action potentials in cardiac and neuronal 
cellsJ^ Another example at the cellular level is the cal- 
cium waves that occur inside cells during critical develop- 
mental events such as fertilization or changes in motility. 
Slower examples include changes in population dynamics 
that occur during disease epidemics and the ignition of 
vegetation in forest fires. Students can use theoretical 
and electronic excitable systems to explore the dynamics 
of these situations and test proposed control mechanisms 
such as cardiac defibrillation, vaccine intervention, and 
fire-fighting strategies. In this paper we model reentrant 
tachycardia responsible for atrial flutter, a common car- 
diac arrhythmia, and its treatment with surgical ablation 
(destruction of the abnormal tissue) . 

One approach is to start with a mathematically simple 



excitable system and then build its electronic circuit real- 
ization. For this approach we use the Fitzhugh-Nagumo 
equations. ^'"^ Another approach is to start with a simple 
circuit that exhibits excitable behavior; we use a three 
transistor excitable circuit. We then find its mathemat- 
ical representation. In both cases the models are pre- 
sented in the form of an excitable system with excitatory 
variable u and inhibitory variable v; 



du 
dv 



(la) 
(lb) 



Excitation pulses are supported by the reverse A^-shape 
of the M-nullcline, f{u, w) = in phase space, and by 
e ^ 1 causing a slow response of v. 

Bunton et al4 describe an excitable circuit based on 
comparators and multiplexer integrated circuits and re- 
view other circuit models. Part of the appeal of their cir- 
cuit is its simplicity, whose behavior can be understood 
by the action of a four-position switch without the use 
of equations. In this paper we emphasize using models 
in the form of Eq. ([T]) to predict the behavior of cir- 
cuits. The level of mathematics is similar to that found 
in current investigations into cardiac arrhythmias i^"— Nu- 
merical solution of the differential equations is easy with 
today's desktop computers using Python or many other 
popular programming environments. The projects de- 
scribed here are intended to introduce the behavior of 
excitable membranes and their actions in cardiac and 
neural systems. The projects also provide valuable ap- 
plications of computational physics and circuit analysis 
and construction. 

Under normal conditions, action potentials that regu- 
late heartbeat travel through the cardiac conduction sys- 
tem in a unidirectional fashion from the sinoatrial node 
to the atrioventricular node, through the atrioventricu- 
lar bundle of fibers and eventually through the purkinje 
fibers in the wall of the ventriclesJ^ This cardiac con- 
duction system can malfunction and lead to a condition 



known as reentrant tachycardia, in which abnormally 
rapid heartbeats circulate through the conduction system 
without stopping. Several types of medically important 
tachycardias exist. The behavioral effects include short- 
ness of breath, dizziness, fainting, heart palpitations, and 
chest pain. Tachycardias may occur in the absence of any 
obvious disease conditions, or may result from an over- 
active thyroid, electrolyte imbalance, or prior myocardial 
infarction (heart attack). Tachycardias may be treated 
with drugs, physical interventions, or in serious condi- 
tions, ablations of the abnormal tissue)^ 

In Sec. II we briefly describe the fundamentals of ex- 
citations in cardiac and neural tissues and then discuss 
the Fitzhugh-Nagumo equations and associated circuit. 
In Sec. Ill we present the three-transistor excitable cir- 
cuit followed by its two-variable model. Section IV gives 
an example of modeling a common cardiac arrhythmia, 
reentrant tachycardia, and its treatment by surgical ab- 
lation. 



II. BIOLOGICAL EXCITABILITY AND THE 
FITZHUGH-NAGUMO MODEL 

In neural and cardiac tissue the momentary excitation 
pulses of the membrane potential (voltage difference from 
interior to exterior across the cell membrane), known as 
action potentials, are caused by ionic currents that flow 
through channels in the cell membrane. Fast-response 
positive feedback arises from the increased conductance 
through sodium and/or calcium ion channels that allow 
those ions to flow into the cell, thereby depolarizing or re- 
ducing the membrane potential. (The resting membrane 
potential is negative, so an inward current depolarizes 
the membrane by raising the potential toward zero.) In 
cardiac cells calcium ion currents cause the cardiac ac- 
tion potentials to have longer durations than most neu- 
ronal action potentials. The slow inhibitory response is 
caused by a delayed opening of potassium channels that 
allow those ions to flow out of the cell, thereby causing 
the recovery or repolarization (hyperpolarization) phase. 
A slowly responding decrease in the conductance of the 
sodium channels also adds to the recovery phase. The 
numbers and types of ion channels vary depending on 
the particular type of cell. The types of ion channels are 
distinguished by their voltage-gated threshold value, ki- 
netics, and ion selectivities. Various ion pumps act like 
batteries to maintain ion concentration gradients across 
the cell membrane. 

The Hodgkin and Huxley model for the neuronal mem- 
brane potential accounts for many of the details of the 
sodium and potassium channel activities,^" and accu- 
rately reproduces the shape of the action potential, at 
the expense of involving four variables. The two- variable 
Fitzhugh-Nagumo model is simpler and has excitable dy- 
namics similar to the Hodgkin and Huxley model. As a 
result, Fitzhugh-Nagumo and similar models are popular 
in studies of cardiac and neuronal dynamics. 
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FIG. 1: Numerical computation (lines) and measured circuit 
data (dashed) for a self-exciting Fitzhugh-Nagumo system 
with source term s{t) = 0.06. The excitatory variable u has 
a large amplitude, and the inhibitory variable v has a small 
amplitude. The voltages have been normalized and the time 
t is dimensionless (see Table [l|. 



The Fitzhugh-Nagumo equations expressed in the form 
of Eq. H]) are^'^ 



du 

dv 
'di 



u{u - a)il - u) - V + s{t) (2a) 

e(w - bv). (2b) 



Typical values of parameters that result in excitable be- 
havior are a — 0.15, e — 0.01, and b — 2.5. The 
threshold-triggered positive feedback is due to du/dt be- 
coming positive for a < u < 1 while v and s{t) are both 
zero. The slow response of the inhibitory variable v is 
ensured by the small value of e. The function s{t) is a 
source term. For example, setting s{t) = 0.06 results 
in self-stimulating periodic excitation pulses as shown in 
Fig. 1. Alternatively, a time dependent source s{t) such 
as a pulse train can be used to stimulate excitation pulses. 

Figure 2 shows the phase space trajectory of the ex- 
citation pulses in Fig. 1 and the nullclines of Eq. ^. 
The M-nullcline is nonlinear and has the reverse- A^ shape 
characteristic of systems that support excitation pulses. 
Fixed points occur at the intersection of the nullclines 
where the time derivatives of u and v are both zero. For 
s{t) — 0.06 the fixed point is unstable resulting in the 
cyclic phase space trajectory. For s{t) = the nullclines 
intersect at the origin, and thus it is easy to use the Ja- 
cobian of Eq. ([5]) in a linearization stability analysis to 
show that the origin is a stable fixed point characterized 
as a spiral nodeJ^ For a proper choice of initial condi- 
tions for s{t) —Q there is one excitation pulse before the 
system settles to the fixed point. 

The concepts of membrane potential and ion currents 
lead naturally to thinking about these excitable tissues 
as electronic circuits, van der Pol and van der Mark sug- 
gested that the heart's dynamics can be modeled with a 




FIG. 2: Numerical computation of phase-space plot (line) 
and nuUclines (it-dashed, w-dot-dash) for the self-exciting 
Fitzhugh-Nagumo system. 



negative-resistance relaxation oscillatorjii Nagumo's cir- 
cuit based on the Fitzhugh-Nagumo model has dynamics 
that are analogous to action potential propagation.'^ Ref- 
erences |ll| and y inspired many publications using elec- 
tronic excitable circuits to investigate cardiac and neural 
systems. For the two circuits we present, a capacitor rep- 
resents a piece of the cell membrane and the capacitor 
voltage is analogous to the membrane potential. Cur- 
rents that charge the capacitor are analogous to inward 
membrane ion currents. Currents that discharge the ca- 
pacitor are analogous to the outward potassium mem- 
brane currents. Voltage thresholds and time constants in 
the circuits account for the ion channels' voltage-gated 
behaviors and kinetics. 

In the Fitzhugh-Nagumo circuit that we use (see 
Figs. 3 and 4) the capacitor voltage V and inductor cur- 
rent / are the variables analogous to u and v in Eq. ^. 
The circuit block V{V) produces voltage V which is a 
cubic function of the voltage V. The input impedance to 
the V{V) block is high so that the current at its input 
may be assumed to be zero. The differential equations 
for V and / are found by applying KirchhofF's current 
law to the capacitor, and the voltage law to the branch 
containing the inductor. 



C 



dV _ V~V 

dts 
dl 



R 



I + L 



dts 



-(V-IRl). 



(3a) 
(3b) 



Here ts indicates that the time variable has units. Equa- 
tion ([3]) is made equivalent to Eq. ^ by using s{t) — 
IsR/2.5, b — Rl/R, and converting to dimensionless 
variables u, v, t, and e defined in Table H] The factor 
of 2.5 is the scaling factor of the MLT04 multiplier IC. 
The voltage V is constructed so that V — V corresponds 
to the nonlinear term u{u — a)(l — u). Note that e is the 
ratio of the time constants RC and L/R for charging the 




FIG. 3: Schematic of the Fitzhugh-Nagumo circuit. R = 
100 rj, C — 0.1 /if, and L — 0.1 H. Rl accounts for the in- 
ductor's intrinsic resistance (220 51), and an external resistor 
(33 0). The function block V{V) is shown in Fig. 4. 



capacitor and establishing the inductor current, respec- 
tively. Because e must be small for excitation pulses to 
occur, we obtain the conditions that charging the capaci- 
tor must be the fast process and creating the discharging 
current in the inductor is the slow one. The analogy to 
cardiac cells and neurons suggests that u{u — a){l ~ u) 
represents the inward ion currents and v represents the 
outward ion currents. 



Dimensionless 


Fitzhugh- 




Quantity 


Nagumo 


Three Transistor Circuit 


u 


1//2.5 


1//5 


V 


iR/2.5 


Vb/5 


t 


ts/RC 


ts/RfC 


e 


R^C/L 


RfC/RACA 



TABLE I: Conversions from measured voltages and currents 
to dimensionless quantities. For the Fitzhugh-Nagumo cir- 
cuit, Fig. 3 defines V, /, R, C, and L. For the three transistor 
circuit. Fig. 7(a) defines V, Vb, Rf, C, Rsi, and C'si- 



Component values are chosen so that e « 0.01 and 
b « 2.5. To avoid the speed limitations of the op 
amps and the MLT04, we choose RC = 10 /zs by using 
R = 100 f7 and C = O.lO^f. The inductor has intrin- 
sic resistance of 2.2ri/mH so an inductor resistance of 
Rl — 250 51 corresponding to 6 = 2.5 gives L = 0.114H. 
We include an external resistance in series with the induc- 
tor so that the inductor current can be measured easily 
from the voltage across this resistor. We use 33 il which 
leaves about 220 57 for a nominal value L = O.IOH; Rl 
in Fig. 3 accounts for both the inductor's resistance and 
the external 33 57. 

The voltage V^ is a function of V. Both voltages need 
to be normalized by the 2.5 scaling factor of the MLT04. 
With u = F/2.5 the required function for u ~ V/2.5 is 

u = u{u — a)(l — u) + u = —u{u — ri){u — 7'2)- (4) 

For a = 0.15 the roots ri and r2 are 1.66 and —0.512. 
Figure 4 shows the circuit used to create the function 




FIG. 4: Circuit for generation of function V{V). Only two of the MLT04's four multipliers are used. The unused multipliers 
(not shown) should be grounded. 



V{V) from Eq. (|3]). The roots are multiplied by 2.5 giv- 
ing the values —4.15V and 1.28V shown Fig. 4. Their 
associated polynomial factors are divided by 2 at the in- 
verting addition amplifiers prior to the multiplier IC to 
bring the voltages within the ±2.5 V range of the MLT04. 
Therefore the result after the multipliers is amplified by 4 
at the final inverting amplifier. We use potentiometers to 
set the values —4.15 and 1.28 V, making it easy to change 
the roots in Eq. (JH). Comte and Marquie use a circuit 
with four fewer op amps to create a cubic polynomial!^ 
This reduction in circuit components comes at the ex- 
pense of not being able to adjust the function so easily. 
They use their circuit to simulate nerve action potential 
propagation^^ 

We used inexpensive op amps (LF412) and multiplier 
IC (MLT04) to demonstrate that high precision compo- 
nents are not required to obtain the results shown in 
Fig. 1. Nominal component values were used for the re- 
sults shown in Fig. 1. The measured capacitor voltage 
shown in Fig. 1 was normalized by the MLT04's scaling 
factor of 2.5. The inductor current was obtained from 
the voltage across the external 33 fl and then converted 
to the dimensionless v by multiplying by 100J7/(2.5 V). 

Nagumo's circuit based on the Fitzhugh-Nagumo equa- 
tions uses a tunnel diode for the nonlinear element;^ We 
use the V{V) block instead, because tunnel diodes are no 
longer inexpensive and readily available. The V{V) block 
has the added benefit of allowing us to easily change the 
functional form of the nonlinear term. 

A constant current source for Is is provided by the 
circuit shown in Fig. 5. The current Is is the voltage 
drop across Re divided by Re- The emitter voltage is 
approximately 0.65 V above the base voltage due to the 
forward biased emitter-base junction. As an example, 
i?i = 20kri and i?2 = lOkil gives a base voltage of 3.3 V 



and emitter voltage of 4.0 V. For Re = 1 kfJ the current 
is 1.0mA. These values correspond to s{t) = IsR/2.5 = 
0.04 in Eq. (Ha} for i? = 100 f^ in Fig. 3. Use of a 25kri 
potentiometer for i?2 allows s{t) to be in the range to 
0.085. 

A pulse train for the current source Is is provided by 
the circuit in Fig. 6 based on a 555 timer. The train's 
period and each pulse's width and amplitude are con- 
trolled by Ra, Rb, and Re, respectively. For example, 
with Ra = 15kn, Rb = 47017, and Re = 680 f7 the pe- 
riod, width, and height are approximately 1.1 nis, 35 /iS, 
and 2.5 mA, respectively. This current pulse corresponds 
to s(t) making transitions between and 0.1. Each pulse 
in this train initiates an excitation pulse in the Fitzhugh- 
Nagumo circuit. A reduction in the train period to 750 /iS 
results in "missing" excitation pulses, a condition known 
as alternans. In cardiac research, alternans induced by 



+5 



-5 



FIG. 5: Constant current source Is- Ri = 20kn, R2 = lOkil, 
and Re = 1 kO, gives Is — 1.0 mA. The transistor is a general 
purpose pnp such as 2N3906. 
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FIG. 6: Pulse train current source Is. Ra ~ 100 kfJ po- 
tentiometer, Rb = 2kf2 potentiometer, C = 0.1 /if, Ri — 
R2 = 22 kn, and Re ~ 1 kf2 potentiometer. The pulse pe- 
riod, width, and amplitude are controlled by Ra, Rb, and 
Re, respectively. 



external pacing provide information about heart func- 
tion. 

Figure 3 looks similar to a RLC circuit. However 
its behavior is not that of a damped harmonic oscil- 
lator. In a harmonic oscillator, the current goes back 
and forth through the capacitor and inductor. In the 
Fitzhugh-Nagumo circuit the capacitor is charged by cur- 
rent through R and discharged by current through L. 



III. THREE TRANSISTOR MODEL 

A simple three-transistor circuit that has excitable be- 
havior similar to the Fitzhugh-Nagumo system is shown 
in Fig. 7(a). The capacitor voltage V and the transistor- 
base voltage Vb are the excitatory and inhibitory vari- 
ables analogous to u and v. This circuit is much simpler 
than the Fitzhugh-Nagumo circuit in Figs. 3 and 4. It is 
possible to calculate all voltages and currents of the cir- 
cuit (using SPICE for example). However, our goal is to 
use the excitable two- variable model in Eq. ([T]). There- 
fore, we model the circuit by the simpler circuit shown 
in Fig. 7(b) in which the two transistors responsible for 
the fast response positive feedback are replaced by the 
resistance i?fast(^) which is a function of the capacitor 
voltage V. The analysis of the circuit in Fig. 7(b) results 
in equations for the voltages V and Vb, 



dV _ 5-V 5-V 

dfs " Ri^stiV) ^ Rs 



V 



lookn 



a 



dtf, 



V~Vb 

Rs\ 
V -Vb 

Rs\ 



-Ic{V,Vb) 

~lB{V,Vb). 



(5a) 
(5b) 



Here, Ic and Ib are the collector and base currents 
for the transistor. This transistor provides the slow in- 
hibitory response and is modeled using the Ebers-Moll 
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(a) Complete circuit 
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(b) 



V 

100k 



V, 



X 



(b)SimpIe model 



FIG. 7: (a) Three-transistor excitable circuit and (b) its sim- 
plified model in which the two transistors responsible for 
the fast-response positive feedback are replaced by -Rfast(V^). 
Rf = IkD., C ^ 0.33 Aif, Ca = l/^f, and Ra = 33k!r2. Rs 
controls the period of self-exciting pulses. The npn and pnp 
transistors are 2N3904 and 2N3906 respectively. 



equationsJ^ii^ The slow response is determined by Rs\ 
and Csi at the transistor's base. The current (5 — V)/ Rg 
is similar to the source Is in Eq. ([5a|l . If Rs is removed 
{Rs — 00), this is analogous to setting Is = 0, and 
s{t) — in Eq. ((2a|) . A wide range of periods for self- 
exciting pulses can be obtained by varying Rs from about 
10 kn to 1 MVL. 

Equation ([5]) is put into the form of Eq. ^ by using 
conversions in Table IH The result is 

du / \ / \Rf Rf 



dv 
'di 



Rsi 5 

Rs\ J , . 
5 



(6a) 
(6b) 



^ — • • i i 




Dimensionless Voltage u 



Dimensionless Time 



FIG. 8: Normalized measured conductance (filled circles) 
of the two-transistor fast-response positive feedback part of 
the three-transistor excitable circuit, and its calculated value 
(line) Qfastiu) from Eq. ([S]). The conductance is normalized 
by 0.001 n~^, the conductance of Rf in Fig. 7(a). 



FIG. 9: Numerical calculation (solid line) and circuit mea- 
surement (dashed line) for the three-transistor excitable cir- 
cuit with Rs = 330 kJ7. The excitatory variable u demon- 
strates threshold-triggered positive feedback, rising slowly 
from to about 0.1 where it then rises rapidly. The inhibitory 
variable v undergoes smaller variations near 0.1. 



where gtastiu) — Rf/RfastiV)- Note that e = 0.01 for 
our choice of circuit parameters, indicating that charging 
the capacitor is the fast process and increasing the tran- 
sistor's base voltage (thereby increasing the discharging 
current) is the slow one. Note that the terms 5fast(u), 
/c(w, w), and Ib{u,v) make Eq. ^ more complicated 
than the Fitzhugh-Nagumo system, Eq. ([2]). 

In the Appendix we show that the collector and base 
currents are 

/p(y, y) = _^ ^e200(.-n) _ 1] + J^ \^^2mv _ g200(.-«)] 

(7a) 
Is{u, v)^^ [e^oo'" - 1] + ^ [e200(-") _ l] . (7b) 

Pf Pr 

The SPICE values from the data sheet for a 2N3904 
are /3/ = 416, ^r = 0.737, and /q = 6.7 x lO^^^A. 
A suitable model for 5fast(w) is found by inspection of 
the fast-response portion of the circuit in Fig. 7(a). As 
V increases from zero and passes through some thresh- 
old value, the two transistors turn on, allowing current 
to flow through the series combination of Rf and the 
pnp transistor. This current provides the threshold- 
triggered positive feedback that rapidly charges the ca- 
pacitor. Thus i?fast(^) needs to have a steep transition 
from very large resistance down to Rf . The inverse of the 
resistance is the conductance, which must make a transi- 
tion from zero to nearly l/Rf. We construct a function 
with these properties and fit it to data. The function we 
use is 



fffast(w) 



R 



f 



1 



RiMV) [1 + e«'i(^t"-5«)] ^1 + (l|ia)-2] 



(8) 

Figure 8 shows Eq. ([5]) with the parameters Wi — 30, 
Vthi = 0.48 V, W2 = 3.5, and T4h2 = 1.25 V, along with 



the experimentally determined conductance obtained by 
measuring the current through Rf (and the pnp transis- 
tor) as a function of the capacitor voltage V = 5u. 

Figure 9 shows measurements of the normalized volt- 
ages u and V with Rs — 330 kfi, along with predictions 
using Eq. ([5]). The phase-space plot in Fig. 10 shows the 
reverse-iV shape for the u-nuUcline that supports exci- 
tation pulses. The nullclines were found by setting the 
derivatives in Eq. ([6]) to zero and numerically finding the 
zeros of v as a function of u. 
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FIG. 10: Numerical calculation of the phase space plot (line) 
and nullclines (u-dash, u-dot-dash) for the three transistor 
circuit. 



IV. APPLICATION TO ATRIAL 
TACHYCARDIA 

As an example of the type of system that can be in- 
vestigated with these circuits, we model atrial flutter, a 
tachycardia (fast heart rate) caused by reentrant circu- 
lation due to unidirectional block, and its treatment via 
surgical ablation. We include descriptions of the relevant 
physiology of the heart and the conditions leading to this 
arrhythmia. 

The self-exciting systems whose behavior is shown in 
Figs. 1 and 9 are analogous to the sinoatrial node, the 
natural pacemaker of the heart, located at the top of the 
right atrium. The source terms s(t) and (1 — u)Rf/Rs in 
Eqs. (|^ and (pa)) correspond to the calcium ion currents 
that flow into cells of the sinoatrial node. These currents 
are responsible for the rising membrane potential which 
is the leading edge of the cardiac pacemaker pulse. The 
heartbeat signal generated by the sinoatrial node prop- 
agates through the contractile muscle cells of the atria, 
causing the atria to contract and send blood into the 
ventricles. Simultaneously, the signal travels through the 
modified, non-contractile cells of the cardiac conduction 
system to arrive at the atrioventricular node located be- 
tween the right atrium and ventricle. Here the signal is 
delayed and ultimately sent down the specialized electri- 
cal conduction cells known as the Bundle of His to ini- 
tiate the ventricular contraction. The delay allows time 
for blood to be emptied from the atria into the ventricles 
before contraction. 

Heartbeat signals (and nerve signals) propagate as ac- 
tion potentials in an excitable medium. This medium is 
modeled by coupling excitable circuit cells together with 
resistor Rd as shown in Fig. 11, similar to Nagumo's cir- 
cuit based on tunnel diodes.'^ Terms must be added to 
the right side of Eq. ([5a|) [or Eq. ((5a)) ] to account for cur- 
rent flow to and from neighboring circuit cells. These two 
terms are rearranged into the form of a diffusion term. 



Vi+i-v, y. -i/.-i 



v^ 



i+l 



2Vi + V,-i Ax'^ d^V 



Rd 



Rd 



Rd 



Rd dx^ ' 
(9) 
where Ax is the distance between centers of neighbor- 




FIG. 11: Three coupled excitable cell circuits that simulate a 
l-dimensional excitable medium. The characteristic time for 
diffusion between neighboring cells is RdC. 



ing cells of the medium. Diffusion is characterized by 
the diffusion coefficient D where the characteristic time 
for diffusing the distance Ax is Ax / D. Thus, the con- 
nection between diffusion and the electrical parameters 
is 



Ax' 
D 



RdC. 



(10) 



The characteristic time for diffusion between neighbor- 
ing cells corresponds to the RC time constant for the 
capacitor and coupling resistor Rd- 

Because multiple excitable cells are required, we choose 
the simpler three transistor circuit in Fig. 7(a). Six of 
these cells are connected in a ring coupled by Rd — 47 kO. 
A self-firing cell with Rg — 330 kfi represents the sinoa- 
trial node. The other cells do not have source resistor 
Rg. The cell on the ring opposite to the sinoatrial cell 
represents the atrioventricular node. The two sides of the 
ring represent two paths from the sinoatrial node through 
the atrial tissue to the atrioventricular node. Excitations 
started at the sinoatrial cell propagate down both sides 
of the ring and meet at the atrioventricular cell where 
they terminate because they cannot propagate through 
each other due to their refractory phases. 

The corresponding predictions provide a challenging 
computational physics project because they involve the 
numerical solution of partial differential equations. We 
use a simple Euler method with a time step At = 2.5 /xs 
chosen to be much smaller than the characteristic times 
for the positive feedback, RfC — (1 kri)(0.33^f) = 
330 Ats and for diffusion, RdC ^ {47 kn){0. 33 ^li) = 
16 ms. Figure 12 shows the numerical predictions and 
data for voltage at a cell along one side of the ring as ex- 





0.04 0.05 

Time (sec) 



FIG. 12: Numerical calculation (solid line) and measurement 
(dashed line) of the capacitor voltage at one cell on one side 
of the 6 cell ring in which a source cell initiates pulses that 
propagate down both sides and meet at the cell opposite to 
source. The period is about 18 ms. This behavior represents 
normal propagation of the heartbeat signal from the sinoatrial 
node through multiple paths merging at the atrioventricular 
node. 
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FIG. 13: Coupling for either normal propagation (SI and S2 
closed), unidirectional block (SI closed, S2 open), or ablation 
(SI open). With SI closed and S2 open pulse propagation 
from left to right is blocked, but from right to left still occurs. 



citations generated at the sinoatrial cell propagate down 
both sides to the atrioventricular cell. The pulses ter- 
minate at the atrioventricular cell. Thus the measured 
period is the same as the firing period of the sinoatrial 
cell, about 18 ms. 

Unidirectional block occurs when propagation is in- 
hibited in the usual (orthograde) direction, but occurs 
for the backward (retrograde) direction. Thus, unidirec- 
tional block in one side of the ring results in a circulating 
pulse because the two pulses no longer meet and cancel 
at the atrioventricular cell. If the circulation period is 
less than the firing period of the sinoatrial cell, then the 
heartbeat rate at the atrioventricular cell will be higher 
than normal, which results in the tachycardia known as 
atrial flutter. In the ring circuit unidirectional block is 
implemented by including switches and a diode in series 
with one of the coupling resistors as shown in Fig. 13. 
Figure 14 shows our numerical predictions and data for 
the voltage at a cell for an excitation circulating around 
the 6 cell ring due to unidirectional block. Note that the 
period is 11ms, less than for the normal signal propaga- 
tion in Fig. 12. In reality, the atrial tissue may contain 
numerous paths of different lengths involved with reen- 
trant dynamics. 

The normal beat period apparent in Fig. 12 is restored 
by breaking the circuit in the path containing the unidi- 
rectional block by opening switch SI in Fig. 13. Breakage 
of the circuit corresponds to ablation (typically using ra- 
diofrequency heating) of the tissue containing the path- 
way with the unidirectional block. 



SUMMARY 



We have presented two excitable systems, and an ex- 
ample of modeling reentrant tachycardia and its treat- 
ment by surgical ablation. These excitable systems may 
be used in a variety of explorations of excitable dynam- 
ics. Suggested possibilities include coupled self-firing ex- 
citable cells, external pacing (using Fig. 6) of an excitable 
media, and determining the minimum path length for an 
excitation circulating on a ring. 
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FIG. 14: Numerical calculation (solid line) and measure- 
ment (dashed line) of the capacitor voltage at one cell on 
a 6 cell ring with circulating excitation caused by unidirec- 
tional block. The circulation period of about 11ms is shorter 
than in Fig. 12. The increased frequency due to circulating 
excitations is manifested as the reentrant tachycardia. 
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Appendix: The Ebers-Moll model for transistors 

The Ebers-Moll model for transistors is often unfamil- 
iar to physics students. The advantage of this model is 
its applicability to all situations for the transistor, from 
fully off with base voltage equal to emitter voltage, all 
the way to fully on (saturation) with the collector volt- 
age nearly equal to emitter voltage. Thus this transis- 
tor model covers both linear amplification and switching 
applications. We give a brief explanation of the model 
assuming familiarity with the basics of bipolar junction 
transistors. 

The standard current-voltage relation for a pn semi- 
conductor junction is 



nv) - /o 



exp 



\Vt) 



1 



(A.l) 



where Vt — kT/e w 25 mV at room temperature and /q 
is a reverse saturation current. For an npn bipolar junc- 
tion transistor the base-emitter and base-collector form 
opposing pn junctions. A voltage across one of the junc- 
tions causes current given by Eq. (|A.ip . The transport 
factor a is the fraction of this current that diffuses across 
the base region thereby flowing through the second junc- 
tion. In addition, the voltage across the second junction 
causes current at that junction given by Eq. (jA.l[) . This 



current has an associated transport current at the first 
junction given by a transport factor for that direction 
of current flow. The result is that the transistor's col- 
lector and emitter currents each have two contributions, 
one similar to Eq. (jA.l[) due to the junction voltage at 
that terminal and the other due to the transport current 
associated with the voltage at the other junction. The 
difference between the collector and emitter currents is 
accounted for by the base current. Ebers and Moll show 
that the collector and emitter currents are-*^^ 
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(A.2a) 
(A.2b) 



Each transport factor is related to the more familiar cur- 
rent gains Pf and /3r by /3 = a/(l — a). A transistor 
is specified by /3/, /3r, and /q. Equation (jA.2p is used 
to express the collector and base currents for the slow- 
response transistor in Fig. 7(b) in terms of the capacitor 
voltage V, base voltage Vb, and the transistor parame- 
ters. The base current is Ib ~ Ie ~ Ic- Equation ([7]) 
is obtained by noting that Vbe ~ Vb, Vbc = Vb — V, 
1/Vt — 40 V^^, and converting to dimensionless u and 
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